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Abstract 

A high-resolution calorimetric spectroscopy study has been performed on pure glycerol and 
colloidal dispersions of an aerosil in glycerol covering a wide range of temperatures from 300 K 
to 380 K, deep in the liquid phase of glycerol. The colloidal glycerol+ aerosil samples with 0.05, 
0.10, and 0.20 mass fraction of aerosil reveal glassy, activated dynamics at temperatures well above 
the T g of the pure glycerol. The onset of glass-like behavior appears to be due to the structural 
frustration imposed by the silica gel on the glycerol liquid. The aerosil gel increases the net viscosity 
of the mixture, placing the sample effectively at a lower temperature thus inducing a glassy state. 
Given the onset of this behavior at relatively low aerosil density (large mean-void length compared 
to the size of a glycerol molecule), this induced glassy behavior is likely due to a collective mode 
of glycerol molecules. The study of frustrated glass-forming systems may be a unique avenue for 
illuminating the physics of glasses. 

PACS numbers: 61.43.Fs, 65.40.Ba, 65.60.+a, 67.40.Fd 
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I. INTRODUCTION 



The glass-forming liquids are interesting materials for study and have continually drawn 
the attention of researchers because of their unusual material properties 
These materials show different thermal transitional and dynamical effects when they enter 
the glassy state. That is, the glassy state is governed by the disorder of the molecules when 
quenched rapidly below the glass transition temperature (T ). Much attention has been 

u n n 

devoted to the study of pure glass-forming liquids such as glycerol |2J, |6|, |]| since the glassy 
state of the pure liquid can be achieved with ease just by cooling below T 9 (^ 190 K for 
glycerol). It remains unclear whether the glass state is due to molecular interactions or 
a collective phenomena of many molecules that prevents crystallization and "freezes" the 
sample into an amorphous state. One avenue to illuminate this phenomena has been the 
study of mixtures of glass-forming liquids or other solvents to isolate any specific molecular 
interactions. However, only a few such studies have been performed to date, and these have 
not settled the central question p. 

An extension of the idea of using mixtures would be to employ a colloidal dispersion in a 
glass-forming liquid to introduce structural frustration. Such a gel system would introduce 
frustration in a controlled way by simply increasing the density of the colloidal particles. 
This would have the advantage of increasing the effective viscosity of the mixture and would 
have an associated mean void length. The perturbations observed for a given mean void 
length could then be compared to a typical molecule size of the glass-forming liquid to 
determine the extent of any collective mode. The colloidal gel could be obtained using 
aerosil particles, a technique that has already found use in the study of quenched random 
disorder at liquid crystal phase transitions and has been shown to involve coupling of the 
aerosil to the dynamics of the host fluid 0, . 

The most common technique for studying different transitions is differential scanning 
calorimetry (DSC) and modulation differential scanning calorimetry (MDSC) techniques, 



which have been widely used in the study of glasses [Hi, |12| • However, the necessary fast 
scan rates severely distort the behavior of slowly relaxing systems. To measure the heat 
capacity directly and over a wide range of frequencies under near equilibrium conditions, 
ac or modulation calorimetry (MC) has proven to be of great benefit for thermodynamic 
studies |ll 0, Id, Q Our interest is to use a MC technique on frustrated glass-forming 
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liquid systems to probe energy relaxations via calorimetric spectroscopy. 

This paper presents the results of a calorimetric spectroscopy study of colloidal disper- 
sions of aerosil particles in glycerol as a function of aerosil density. The addition of aerosil 
introduces glass like dynamics with relaxations on the order of a second beginning at void 
length scales on the order of 100 nm. This indicates that the induced glassy behavior is 
due to a collective mode of many glycerol molecules. A theoretical review of calorimetric 
spectroscopy and the relevant glass dynamics is presented in Sec. |Hj and the experimental 
details are described in Sec. IIHI The results for temperature and frequency scans as well 
as an Arrhenius analysis are shown in Sec. IIV[ with discussion and conclusions drawn in 
Sec. El 

II. THEORETICAL DESCRIPTION 
A. Complex Heat Capacity 

Since the source of any sample's heat capacity are the fluctuations of the sample's energy, 
it is natural to consider C as a dynamic quantity. However, for most materials the charac- 
teristic relaxation time of the energy fluctuations are too short to be sensed by traditional 
calorimetric techniques. Nevertheless, a complex heat capacity may be defined containing 
a real (C) and imaginary (C") component indicating the storage (capacitance) and loss 
(dispersion) of the energy in the sample in a direct analog to a complex permittivity. Us- 
ing a linear response approach, the relaxing enthalpy fluctuations (at constant pressure) 
denoted as 5Hr defines an enthalpy correlation function as (SHn(0)6Hit(t)). The complex 
heat capacity at a given frequency ui is then given by 

C p (uj) = C' p (uj) - i / uk B {3 2 (5H R (0)6H R (t))dt (1) 
Jo 

where (3 = (fceT) -1 and the time integral multiplying the imaginary number is C" v . For 
energy dynamics characterized by a single mode of frequency u m , the real part has two 
asymptotic limits, fast (w -> oo or > u m , denoted by C^°) and static (w or < u> m , 
denoted by C° = + kB(3 2 {5H R )) compared to this mode. The imaginary part would 
exhibit a peak at uj m that would be the inflection point of the smooth roll over between the 
two frequency limits of C' p . 
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B. Modulation Calorimetry 



Modulation calorimetry allows one to make frequency dependent C v measurements and 
so access calorimetric spectroscopy. In modulation calorimetry, a small oscillating heating 
power P{t) = P exp(iuit) is applied to the sample+cell and induces an rms temperature 
rise as well as small temperature oscillations. Note, in what follows the cell is considered to 
consist of the actual sample holder and the attached heater and thermometer. For a sample 
of finite thermal conductivity and a sample+cell experimental arrangement of a thickness 
less than the thermal diffusion length, the heat flow equations may be set up considering the 
sample+cell as a single thermal mass (i.e., geometry independent) 16]. The temperature 
oscillation amplitude T ac is given to second order by 

where P is the amplitude of the applied heating power, u is the angular frequency, C = 
C s + C c is the total heat capacity of the sample and cell, and T ac is the amplitude of the 
temperature oscillations. In addition, there are two important thermal relaxation times, 
the external relaxation time constant r e = R e C and the internal relaxation time constant 
T u = T s + T ci "which is the sum of square internal time constants of the sample and cell 
(r s = R S C S and r c = R C C C ). Here, R s is sample's thermal resistance and R e is the external 
thermal resistance to the bath. There is also a phase shift between the applied heat and 
resulting temperature oscillations $ but it is more convenient to define a reduced phase shift 
= $ + 7r/2 since for frequencies below 1/ra, $ ~ —n/2. The reduced phase shift, to the 
same accuracy as Eq. (j2J), is given by 

tan(0) = (cUTe) -1 — LUTi (3) 

where 7* = t s + r c . The different internal relaxation times that enter in Eq. (J2J) and (jSJ) are 
related by t\ = rf + 2r s r c and must be taken into account in order to extract the complex 
heat capacity. 

Noting that 1 + tan 2 (0) = cos~ 2 (0), Eq. (JHJ) can be substituted into Eq. (J2J to give 

1/2 

(4) 



C = C* 



2t c t s lu 2 + 2— + 



COS 2 (0) ' T e 3R e 

where C* = Po/(uT ac ) is defined for ease of notation. Factoring out a cos(0) gives 

C = C*cos(0)/(cj) (5) 



where the total heat capacity is identified as the real part of the complex heat capacity. The 
function f(oo) is then given by 

-1/2 

(6) 

and may be regarded as a correction for a sample and cell having comparable finite thermal 
relaxation times. 

For most materials away from any phase transition, the imaginary part of the heat ca- 
pacity is essentially zero at typical heating frequencies (< 15 mHz). This quantity may be 
derived from Eq. Q by substituting the definitions of the relaxation times and Eq. (JHJ) for 
the heat capacity. This gives 

C" = C* sm{<f))g{uj) - {uRe)- 1 « 0. (7) 

where a frequency dependent correction function is introduced as 

M = 'w( l + £®) < 8 > 

that also accounts for the comparable relaxation times of the sample and cell. 

Typically, the cell's geometry and mass may be controlled and so its thermal relaxation 
time can be made much less than that of the sample, r c < r s . In this case, the two internal 
relaxation times are approximately equal. Taking Ta = 7$, the correction function g(u) does 
not change but the uj 2 term vanishes in f(oo). The frequency dependence for both correction 
functions enter through the reduced phase shift. At sufficiently low heating frequencies 
uj < r" 1 , cf) approaches tt/2 (<& approaches 0) and f(oo) ~ g{w) ~ lj which reduces the real 
and imaginary heat capacities in Eqs. (jHJ) and (j2J) to 

C = C* cos(0) (9) 

and 

C" = C*sin(0) - (u}R e y l . (10) 

where any energy dynamics or dispersion in the sample will be indicated by a non-zero value 
of C". 
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C. Brief Review of Glass Dynamics 



Glass-forming or super-cooled liquids exhibit dynamics due to inhibited structural re- 
laxations that diverge upon approaching the glass-transition until the system falls from 
ergodicity. Glass-formers are characterized by a rapid increase in viscosity (from 10~ 2 poise 
above to 10 15 poise below the glass transition temperature) with decreasing temperature, 
which reflects the rapidly shifting timescale of relaxation dynamics of the liquid. The slow 
molecular motion, known as a-relaxation, has been well characterized on timescales between 
100 ps and longer 17, is| and this stops completely below the glass transition temperature. 
The fast motion (/3-relaxation) is predicted to exist at timescales between 1 and 100 ps [19], 
and persists even in the glassy state. However, /3-relaxation is difficult to determine unam- 
biguously by experiment and its physical meaning is debated. Qualitatively, a-relaxation 
can be thought of as the collective motion of many particles while /5-relaxation is the motion 
of a single molecule rattling within a cage of nearest neighbors. Fundamentally, the issue 
of glasses centers on whether the glassy behavior is due to "stuck states" of a collection of 
molecules or due to molecular frustration preventing the crystallization of the material. 

Quantitatively, glassy dynamics typically obey Vogel-Fulcher-Tammann (or Arrhenius) 
behavior indicating that the dynamics are energetically activated I20I . l^ . According to 
Arrhenius behavior, the relaxation time of glass forming liquids can be given by 

r = r exp(/3A£) (11) 

where r is the relaxation time of some structural fluctuation, tq is the high-temperature limit 
of this fluctuation, AE is the activation energy, and (3 = (/c^T) -1 as usual. This activated 
relaxation mode can become very long and complicate the interpretation of thermodynamic 



measurements. Glycerol has been a particularly we 
liquids that exhibit all of the above characteristics 



1 studied example of such glass-forming 
1. 



III. EXPERIMENTAL 



The pure glycerol obtained from Aldrich was used after carefully degassing at 323 K for 
> 2 hours because of its hygroscopic nature. The pure glycerol has a molecular weight 
of M w = 92.09 g mol -1 , a density of p g = 1.26 g cm -3 , and a nominal glass transition 
temperature of T g ~ 195 K. The hydrophilic type-300 aerosil silica nano particles obtained 
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from Degussa |22| were thoroughly dried at ~ 573 K under vacuum for ~ 2 hours prior to 
use. The specific surface area of the aerosils measured by the manufacturer via BET nitrogen 
isotherms is a = 300 m 2 g" 1 and each aerosil sphere is roughly 7 nm in diameter. However, 
SAXS studies have shown that the basic aerosil unit consists of a few of these spheres fused 
together during the manufacturing process Typically, the hydrophilic nature of the 
aerosils arises from the hydroxyl groups covering the surface and allows the aerosil particles to 
hydrogen bond to each other. This type of bonding is weak and can be broken and reformed, 
which leads to the thixotropic nature of gels formed by aerosils. The aerosil gelation in an 
organic solvent occurs via a diffusion-limited aggregation process resulting in a fractal gel of 
df ~ 2.4 |23j]. See Fig.^ However, pure glycerol is a hydrogen-bonding liquid that may alter 
the dispersion of aerosil from that shown in Fig. ^ Each glycerol +sil colloidal dispersion 
sample was created by mixing appropriate quantities of aerosil and glycerol together with 
spectroscopic grade (low-water content) acetone that was subsequently evaporated slowly 
away. The resulting mixtures were then annealed and degassed at 323 K for 1 hour. The 
resulting samples appear by visual inspection to be very similar to those made in an organic 
solvent like a liquid crystal. In addition to the pure glycerol sample, three dispersion samples 
were prepared with a mass fraction of aerosil of 0.05, 0.10, and 0.20. In addition, these 
samples may be characterized by the conjugate density p$ defined as the mass of aerosil per 
open (glycerol) volume that allows one to determine the mean-void length l = 2/aps |l0j. 
The characterization of the aerosil by ps and lo does not depend on the details of the 
dispersion structure, only that it be random. It is likely that the aerosil dispersion remains 
random in glycerol. See Table |J for a summary of these parameters. 

High-resolution ac calorimetry was performed using a home-built calorimeter at WPI. The 
sample cell consisted of a silver crimp-sealed cup+lid where cup dimensions are ~ 12 mm 
diameter and ~ 0.5 mm thick and lid dimensions are ~ 12 mm diameter and ~ 0.1 mm 
thick (closely matching the dimensions of the heater). The average mass of the silver cell 
was 0.135 g and the sample (glycerol-l- aerosil) was 36 mg and the total mass of cell+sample 
did not deviate by more than ~ 35 mg between the different samples. After the sample was 
introduced into a cell and the lid crimp sealed to the cup under pressure a 120-f2 strain-gauge 
heater and 1-MT2 carbon-flake thermistor were attached. The cell was then mounted in the 
calorimeter, the details of which have been described elsewhere 24|. In the ac-mode, power 
is input to the cell as Poe luJt resulting in temperature oscillations with amplitude T ac and a 
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relative phase shift of 0. From these parameters and the all necessary finite conductivity 
corrections, the real and imaginary heat capacity were determined. All cells studied in this 
work closely matched each other in dimension and mass to better than 5% in order to better 
isolate effects introduced by the aerosil. 

Three experiments were performed: temperature, time, and frequency scans with the 
following protocols. For temperature scans, measurements were made on a freshly mounted 
sample at a constant heating frequency of 15 mHz as the cell first heated then immediately 
cooled between 300 and 380 K using a constant rate of ±2 K fir -1 . After the cell had cooled 
back to the initial temperature of 300 K, a time scan was preformed to monitor the heat 
capacity as it relaxed back to the initial value. Frequency scans were then performed at 
fixed temperatures from 300 to 380 K in steps of 20 K over the frequency range from 1 to 
2000 mHz (i.e., u from 0.0063 to 12.6 s" 1 ). 

IV. RESULTS 

A. Temperature and time scans 

Heating and cooling scans for all samples are shown in the Fig. El From the C p data, 
clear indications of hysteresis are observed for the pure glycerol in that the values of C p on 
cooling do not reproduce the values observed on the heating scan of a freshly loaded sample. 
For the glycerol+sil samples, this hysteresis decreases with increasing aerosil content. This 
effect may be quantified by the area of this hysteresis loop (AHh yst ) that decreases in a 
near linear fashion with the mass fraction of aerosil. This behavior is in marked contrast to 
that of the reduced phase shift. For pure glycerol reproduces itself between heating and 
cooling, but begins to take on increasing hysteresis like that described for C p but with 
increasing aerosil content. 

The measurements made as a function of time, not shown here, upon returning to the 
starting temperature clearly indicate a very slow, nearly linear, relaxation back to the initial 
C p values. For the pure glycerol, this was observed to take 2 to 3 days; while for the 0.10 
dispersion sample, 5 to 6 days were required for the sample to recover. This time increased 
with increasing aerosil content. Clearly, long relaxations persist at these high temperatures 
for glycerol and the character of this relaxation changes with the addition of the aerosil gel. 
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These extremely slow relaxations (on the order of days) are of unclear origin. 
B. Frequency scans 

To probe the slow dynamics of the pure and dispersion samples, frequency scans were 
performed from 2000 mHz down to 1 mHz. The pure glycerol frequency scans exhibited 
a dispersion peak in the imaginary heat capacity at a frequency (u p ) coincident with the 
inflection point of the real heat capacity roll-off. This dispersion peak remained stationary 
at Up ~ 0.5 s -1 as the temperature progressively increased from 300 K to 380 K. In contrast, 
the dispersion peak for the glycerol+sil samples clearly indicate a shifting dispersion peak 
with increasing temperature. An example of the resulting real and imaginary heat capacities 
after all calibrations and internal corrections were made is shown in Fig.|3]for the 0.10 sample 
at five temperatures. The imaginary heat capacity shows a dispersion peak at a frequency 
that again coincides with the inflection point of the real heat capacity roll-off. However, as 
the temperature increases progressively to 380 K, the dispersion peak shifts towards higher 
frequencies consistent with results observed in the glassy state of pure glycerol jlj]. 

The evolution of the C'^/u peak at 300 K as a function of frequency is given in Fig. 0]for 
several aerosil concentrations. The dispersion peak increases in height and shifts to lower 
frequencies with increasing amounts of aerosil. The integration of C'^/uj over the range of 
frequency covered yields the total dispersion heat capacity of this relaxation mode (AC"), 
which linearly increases with aerosil concentration. 

The frequency of the imaginary heat capacity dispersion peak uj p determines a char- 
acteristic relaxation time for this dynamic mode as r r = l/u p . To establish whether these 
relaxations are energetically activated, a semi- log plot of r r against 1000/T was made, which 
would reveal a linear region whose slope is directly related to the activation energy AE. See 
Eq. flllj) . Such a plot is shown in Fig. 03 The pure glycerol does not exhibit activated dy- 
namics and the relaxations for the glycerol+sil dispersions become progressively activated 
at progressively lower temperatures representing the onset of Vogel-Fulcher or Arrhenius 
behavior. It is clear that as the aerosil content increases, the activation energy (high- 
temperature relaxation time) of the glycerol+sil system continuously increases (decreases). 
Table summarizes the results of the dispersion peak integration as well as the Arrhenius 
analysis. 
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V. DISCUSSION AND CONCLUSIONS 



Since the sample and cell configuration was constructed to be nearly identical for all 
the samples studied, the differences observed with the introduction of aerosil is particularly 
significant. This is especially true since the thermal conductivity of the silica is greater 
than the glycerol and its increasing content should have improved the internal thermal 
conductivity of the sample. Also, since both the aerosil and glycerol can hydrogen-bond to 
themselves and each other, the aerosil surface is likely coated by a bound layer of glycerol. 
This would strongly couple the two components of this dispersion. At room temperature and 
above, pure glycerol behaves as an isotropic liquid but still retains some strong molecular 
interactions to induce the hysteresis shown in Fig. |21 The mass of the cell+sample was 
checked repeated at various stages during the scan protocol and there were no significant 
changes. Thus, desorbing and absorbing moisture during the thermal scans cannot account 
for the observed hysteresis behavior. However, this hysteresis relaxes very slowly, on the 
order of days, and so it does not seem to be connected to the dynamics observed by the 
frequency scans (on the order of a second) for the pure glycerol or the glycerol+sil samples. 

As the silica density increases, the effective viscosity of the colloidal mixture increases and 
there is also an increase in the activation energy. Since these experiments were conducted well 
above the pure glycerol glass transition temperature of 195 K, the induced glassy behavior is 
likely due to the aerosil-imposed frustration on a large collection of glycerol molecules. This 
is evident considering that the onset of this effect begins even for the 0.05 sample, which 
has a mean void length scale of ~ 100 nm. The estimated size of a glycerol molecule from 
the specific volume is ~ 1 nm. A summary of the thermodynamic results are given in Fig. El 
and Table |lj Clearly there is a strong correlation of the activation energy AE and the total 
heat capacity loss AC' , both increasing linearly with aerosil content. 

The energy dynamics of glycerol containing a colloidal dispersion of aerosil silica nano 
particles has been studied by high-resolution modulation calorimetry. This has revealed 
a temperature and silica density dependent complex heat capacity. It has been observed 
that under certain experimental conditions of temperature and frequency, this system shows 
interesting frequency dynamics and displays dispersion peaks that shift towards higher fre- 
quency with increasing temperature. These results are consistent with those found for the 
glassy state of glycerol [l|. As the silica density is increased in the isotropic glycerol, larger 
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dispersion peaks are found at lower frequencies indicating an increase in the activation en- 
ergy of this relaxation. This phenomena may be explained in terms of the formation of 
"stuck states" of glycerol in the voids between glycerol-coated silica particles also manifest- 
ing itself in the increase in the sample's viscosity. This work highlights the potential of using 
frustrated glass-forming liquids via colloidal gel dispersions as a means of gaining insight 
into the glass state. Dielectric spectroscopy studies of such systems as well as x-ray analysis 
of the structure of the aerosil dispersion would be of great interest. 
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FIG. 1: Cartoon depicting the fractal gel formed by the diffusion- limited aggregation (hydrogen- 
bonding) of aerosil nano-particles in an organic solvent. However, in glycerol H-bonding between 
the glycerol and aerosil is also possible and may change the "pearl-necklace" nature of the aerosil 
dispersion. In any case, the distribution of aerosil should remain relatively random. 
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FIG. 2: Heat capacity (a) and phase shift (b) observed in temperature scans for glycerol+sil 
mixtures containing 0, 0.05, 0.10, and 0.20 % aerosil fraction by mass. Arrows indicate the direction 
of the scan, all at a rate of 2 K hr _1 , with the cooling scan immediately following the heating. 
All scans were performed using a heating frequency of 15 mHz. Note the hysteresis in C p that 
progressively decreases with increasing silica content. 



15 




2.5 3.0 3.5 4.0 4.5 

co ( rad s~' ) 



FIG. 3: Real (a) and imaginary (b) heat capacity as a function of heating frequency for the 0.10 
sample upon heating a fresh sample to selected temperatures from 300 to 380 K. The lines are 
guides to the eye and inset denotes the different temperatures in kelvin. 
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FIG. 4: Dispersion peak scaled by the frequency for all samples at 300 K. The integration of this 
peak yields the total dispersion or loss heat capacity, ACp. The inset lists the mass fraction of 
aerosil for each sample. 
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FIG. 5: Semi-log plot of relaxation time r = l/u p (taken from the frequency of the peak in Cp) 
versus inverse temperature for each sample studied. The inset lists the mass fraction of aerosil. 
The slope provides the activation energy AE and the intercept at 1000/T = provides the high- 
temperature relaxation time tq of the samples. 
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FIG. 6: Density dependence of the energy dynamics. Panel (a) shows the hysteresis energy 
(AHhyst) on left and activation energy (AE 1 ) on right, and panel (b) depicts the total dispersion 
heat capacity (AC") of samples containing 0, 5, 10, and 20% silica by mass. 
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TABLE I: Summary of the calorimetric results for the glycerol+aerosil samples. Shown are the 
silica mass fraction x = M si i/MT tal, conjugate density ps (in g of aerosil per cm 3 of glycerol), mean 
void length Iq (in nm), integrated dispersion peak AC" (in J K -1 mole -1 ), the high temperature 
the activation energy AE (in kJ mole -1 ), and the high-temperature relaxation time To (in seconds) 
obtained by an Arrhenius analysis. 
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